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Abstract 

The oscillations of pseudo-Dirac neutrinos in matter are discussed and applied to the 
solar neutrino problem. Several scenarios such as both u e and being pseudo-Dirac and 
only v e or being pseudo-Dirac are examined. It is shown that the allowed region in the 
mass-mixing angle parameter space obtained by comparing the solar neutrino data with 
the calculations based on the standard solar model and the MSW effect is not unique. 
The results depend on the nature of neutrinos; for example, if both v e and are pseudo- 
Dirac, the allowed region determined by the current solar neutrino data does not overlap 
with that obtained in the usual case of pure Dirac or Majorana neutrinos. 
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The observation of the solar neutrinos, combined with the standard solar model [lfl 
and calculations of the MSW effect can provide important clues to understanding the 
basic properties of neutrinos. Current data from Homestake ||, Kamioka |4[] and SAGE 
U have already narrowed down considerably the allowed region in the mass-mixing angle 
parameter space. 

Recently there has been a proposal || to explain the solar neutrino problem by using 
the MSW effect with only one generation of pseudo-Dirac |7|] electron neutrinos with a 
large transition magnetic moment. Majorana neutrinos emerge naturally in most exten- 
sions of the standard model, but some models (see, for example, Ref.H) yield pseudo-Dirac 
neutrinos, which were used, among others, to explain the solar neutrino puzzle. 

In this paper, we generalize the one generation picture to the case of two generations 
and examine the consequences. In particular, it is shown that in the pseudo-Dirac neutrino 
case the current data yield entirely different (non-overlapping) allowed regions in the mass- 
mixing angle parameter space from those in the standard Dirac or Majorana neutrino 
cases. 

The Dirac-Majorana mass matrix for the neutrino states v e L and v R is given by 
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(1) 



where Ml and m e L R are Dirac and Majorana masses, respectively. The assumption that 
m e L + m e a = leads to pseudo-Dirac neutrinos generated by a mechanism similar to that 
originally discussed by Wolfenstein 0. In the following we make the assumption that 
Ml 3> \m e L ^ R \ which also leads to pseudo-Dirac neutrinos ||, i.e. two almost degenerate 
(in mass) left-handed neutrino states v\ and i/f (with masses m\ ~ m\ ~ Ml) which are 
expressed as 



v\ = i cos 9 e v e L — i sin 9 e u R 



i>2 = sm8 e i>l + COS^g^Jj 



(2) 



where the factor i guarantees the positivity of the mass eigenvalues. The mixing angle 9 e 
is given by 

2M e 

tan(2tf e ) = - 5 (3) 
(m% - m e L ) 

In our case 9 e ~ 45° because Ml 3> m e L R . Therefore, our pseudo-Dirac neutrinos are 
special ones, with almost 45° mixing angle. In general, however, pseudo-Dirac neutrinos 
can have any mixing angle and whenever pseudo-Dirac neutrinos are generated with one 
generation one has to introduce sterile neutrinos. 
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In the following, we assume that 10 11 eV 2 ^ Am 2 = (m^) 2 — (ml) 2 < 10 7 eV 2 



where the upper limit comes from a cosmological constraint |10[ on the oscillation into 
sterile neutrinos and the lower limit is necessary in order to have a vacuum oscillation 
length much shorter than the sun-earth distance. In this case, during the sun-earth 
propagation, one half of the initial flux of v e L will be depleted due to the maximal (45° 
mixing) oscillations between v e L and ~v% when the time average is taken. Therefore, the 
ratio of the v\ flux at the earth and the initial flux is one half for the SAGE (S), Homestake 
(7i) and Kamioka (/C) experiments: 

S = H = 1C= 1 -. (4) 
If neutrinos have magnetic moments large enough to induce a spin-flip during their 



propagation in the magnetic field of the sun [O, the ratios of the neutrino fluxes are 



e — e — 1 ~2a 1 -2a . . 

K ■ vi ■ < - K = « : ■ — ^— ■ « (5) 

with < a < 0.5. Since a is related to the spin-flip probability P s f as a = (1 — P s f)/2, any 
deviation from a = 0.5 is an indication that spin-flips actually took place. The detection 
rates are then 

1 — 2a 

S = H = a , K = a + 0A2 — - — (6) 

where we have used a{yle~) ~ 0.42cr(z/^e~). Assuming that the standard solar model U 
gives the correct v e L flux produced in the core of the sun, the ratios of the observed fluxes 
and the initial flux are H cxp = 0.27 ±0.04 0, /C cxp = 0.46 ±0.08 g and <S cxp = 0.15 ±0.27 
||. The range of the parameter a for which Eq.(|6]) explains the observed ratios is 0.29 ^ 
a < 0.31. 

Now we generalize the above one generation picture to the case of two generations 
and study its consequences. We assume that, before the mixing between the electron and 
muon sectors, the muon neutrinos are also pseudo-Dirac particles, i.e. ^> so 
that Ofj, ~ 45°. The two almost degenerate (mf ~ m% ~ M£) mass eigenstate muon- 
neutrino states z/f and v% are given by 

Ui = i cos QpV^ — i sin Q^Vr 

_ (7) 
u% = sin 9^ + cos^i/^ 

In the muon neutrino sector, we assume that 10~ n eV 2 < Am 2 = (m^) 2 — (m^) 2 < 
10 _2 eV 2 , where the upper limit is due to the cosmological constraint ||10||. However, both 



3 



the mixing angles among electron and muon neutrinos and Am 2 ~ (Mj 



■ft\2 



(M e D f are 



left unknown. For simplicity, we assume that Am^ ^> Am^ , Am 2 u and that the mass 



eigenstates are given by 
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(8) 



where c e = cos9 efl and s g = sm9 efl . The mixing angle 6 e ^ is practically equivalent to the 
usual mixing angle between v e and in the Dirac or Majorana cases. Also, in Eq.(|8|) the 
pseudo-Dirac mixing angles 6 e and 9^ have been approximated to 45°. The mass matrix 
in the weak basis is given by 



V 



Am 2 2 c 2 e + Am 2 u s 2 
(m| 4 - m\ 2 ) c g s 
{Am 2 u - Am 2 2 ) c g s e 



Am 2 2 c 2 + Am 2 d s 2 



31'" 



{m^ 4 — m 



12 J '-e 



CaSr, 



m 12 c e + m 34 s e 



(Am 



| 4 AlBjj) CgSg 



(AmJ 4 - Amf 2 ) c„s e m\ 2 s 2 + m| 4 Cg + A p 



{m 2 A - m\ 2 ) c e s e 



Am 2 2 s 2 + Am 2 d c 2 



31' 



(Am§ 4 - Am\ 2 ) c e s g 

(7M34 ^12) Ce-^e 

Ara? 2 s* + Am^cJ 



m 12 s e + m 34 c e 



(9) 



with 



m? 2 



Am? 2 



+ m 2 



m 34 = 



m| + m 4 



Arrtg 4 = 



m 



(10) 



A, 



A-cc + A 



NC 



A„, = A 



NC 



A cc = 2V2G F EN e , A NC = -V2G F EN n 



where mi,m2,m3 and m 4 are the mass eigenvalues. The values of the effective mass 
squared in matter are shown in Fig.l for u%, i>2, ^3 and z/ 4 and their antiparticles as 
functions of the matter density p. In Fig.l there are two MSW resonance regions Rl and 
R2 and two possible spin-flip resonance regions Rl m and R2 m in addition to the region 
R discussed in Ref. 0. In the region R the maximal vacuum oscillations lead to the 1/2 
suppression of the v e L flux, as discussed above. 

First we discuss the case in which neutrinos have large enough magnetic moments to 
induce resonant spin-flips and the resonance regions Rl m and R2 m are in the convection 
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zone. In this case the right-handed neutrinos v^i v %i u l and can be generated from 
the original v e L by the resonant spin-flip processes directly or indirectly, e.g. via v e L — + u£, 
u e — > v e — > v e v e — > — > v e 

L L 1 h " Rl L U R L ' 

Regardless of whether the resonance spin- flip processes are adiabatic or not, the flux 
ratios at the earth are expressed as 

l-2a I -2a , 1-2/3 , 1-2/? 

aa : a : a : aa : bp : b : b : bp ( 11 ) 

2 2 H 2 2 H K J 

with a + b = 1 and < a, f3 < 0.5. The detection rates are then 

S = H = aa, /C~aa + 0.21a(l-2a) + ^(l-a) (12) 

where we have used a{v£e~) ~ lcr(z/^e~), and cr(^e _ ) ~ l<j(z/^e~) instead of cx^e") ~ 
l<j(^e _ ) for simplicity. This approximation makes the second equation in Eq.(|12|) free of 



the parameter (3. The results in Eq.([i2D are consistent with iS exp , T~tex P and /C exp for the 
parameter ranges a ^ 0.93 and 0.29 ^$ a ^ 0.31 within la. The fact that a must be very 
close to unity implies that the Landau-Zener transitions at the resonances Rl and R2 are 
extremely non-adiabatic in this model. Since the deviation from a = 0.5 indicates the 
presence of spin-flip, the above range of a requires a large transition magnetic moment 
between v e L and v e R (e.g. ~ lO -10 /^ for B ~ 10 KG in the convection zone). Although there 
exist many models that can yield such a large magnetic moment, they appear somewhat 
unnatural and thus we do not consider this scenario further. 

In the absence of magnetic moments, no right-handed neutrinos are produced as v e L 
(~ z/4 in the core) propagates outward from the core. At the resonance R2, z/4 is split into 

and z/3 with fractions (1 — P R2 ) and P R2 , respectively, where P R2 is the Landau-Zener 
transition probability at the resonance R2. At the resonance Rl, the fraction P R2 of u 3 
is further split into v 3 and z/ 2 with fractions (1 — P R1 )P R2 and P R1 P R2 , respectively, P R1 
being the Landau-Zener transition probability at the resonance Rl. 

Let us assume that the two resonance regions do not overlap so that the two resonances 
can be treated separately. Then one can estimate the respective Landau-Zener transition 
probabilities as follows. 

(1) First resonance (Rl). This is a resonance between v e L and v% and occurs when the 
first and third diagonal elements of are equal, i.e. for A cc = (m^ — m\ 2 ) c 20 . In the 
neighborhood of the resonance the MS W evolution equation is dominated by the v e L -v^ 
2x2 sector: 

m 2 l2 c 2 g + m 2 34 s 2 + A e (m 2 34 - m 2 2 ) 



(m| 4 - m\ 2 ) c e s e m\ 2 s 2 g + m 2 u c 2 + A p 



(13) 
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The Landau-Zener transition probability at the resonance is given by 



cxp 



4 m 2 4 -m 2 2 



(14) 



with h m = i |e 
usual Am 2 . 



The mass factor in the exponent, (m| 4 — mf 2 ), is equivalent to the 



(2) Second resonance (R2). This is a resonance between v e L and ^ and occurs when 



the first and fourth diagonal elements of M 2 are equal, i.e. for A e 



.31 



127 u 29- 



In the neighborhood of the resonance the MSW evolution equation is dominated by the 

L J R 



vf-Up 2x2 sector: 



^ m 12 C e + m 34 S e + A e (AW34 ~ Amf 2 ) C e S e ^ 



M,1 



(Am 2 4 — Am? 2 ) c e s e m? 2 s 2 + m^c^ 



2 „2 



(15) 



/ 



The Landau-Zener transition probability at the resonance is given by 



P R2 = exp 



7T 4 77*4 



12 



Am 2 4 — Am 2 



12 



2 2 

m 34 — m 12 



(16) 



with h. 



l Op 

p cte 



. The ratio of the mass factors in the exponents of Eqs.flU]) and (TfB| ) 



is [(Am 2 4 - Am 2 2 )/(m 2 4 - m 2 2 )] 2 ~ (Am 2 4 /m 2 4 ) 2 ~ (M>^ fl /M^ 



t 2x2 



M \2 



The numerical value of this ratio is supposed to be small for the pseudo-Dirac neutrinos 
under consideration. Therefore, we take the non-adiabatic approximation for the resonant 
transition at R2. In order to see the region of validity for this approximation, let us 
consider the exponent of Eq. (|i~6|) which is written as 



Q 



772 



7rsin 2 (2^) R G Am 



>77 r l 
e\x \ L,R 



4 cos(20 e/ J 10.45 
sin 2 (20 e 



'efi 

-2.6 x 10 3 



E 
Am 2 M 



cos (29, 



(17) 



where we have used E = 10 MeV and for definiteness R /M£ ~ 0.01 which corresponds 
to the mixing angle 9^ = 44.86°. The non-adiabatic region which satisfies {Qi&l < 1 is 
below the solid line in the upper right-hand corner in Fig.2. It will be shown that the 
solution of the solar neutrino problem in the pseudo-Dirac neutrino model indeed lies in 
this region. 
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Therefore, one has the following ratios of fluxes at the earth 



< ■ < 



leading to 



P m cos 2 6 eil + (l-P R1 )sm 2 e, 



eft 



s = n 



[1 - 2P m ) cos(20, 



e/i. j 



K 



P R1 sm 2 9 eil + (l-P R1 )cos 2 9, 



e/j, 



(18) 



1 

12* 



5 

12 



'1 - 2P R 



where we have used the fact that, because of the maximal mixings between v\ 



cos(29 efl 

(19) 



v e R and 



i/r, only half of the neutrinos can be detected. 



Treating P m as a constant parameter (neglecting the energy dependence of the tran- 
sition probability), S, 

values of P m , but are reproduced within 2a for 



ex P , T^exp and /C exp are not reproduced by Eq.(|T^) within la for any 



0.52 < 



~-~(l-2P m )cos(20 e/i ) 



< 0.70 



(20) 



The two-generation pseudo-Dirac neutrinos discussed here are a special case of two- 
generation sterile neutrinos in which 9 e and 9^ are 45°. Since the only mixing angle which 
is relevant in the analysis is 9 e/1 , the above result should be compared with the usual MSW 
effect with two generations of neutrinos with the same mixing angle. In this case one has 



2P LZ ) cos(2#, 



K = I + i 

6 6 



i-Id 

2 2 V 



2P LZ ) cos(20. 



eji ) 



(21) 



with the Landau-Zener transition probability P LZ . Equation ( pT|) can reproduce S exp , TC, 

I_ I(l_2P LZ )cos(2^ 



ex]) 



and /C exp within la if 0.26 ^5 



iS 0.31, and within 2a if 



0.19 < 



\~ i(l-2P LZ )cos(20 eM ) 



< 0.35 



:22i 



There is an important difference between Eqs.fllTiD and fl21~|) : In Eq . (|19D , there is an 
additional factor 1/2 due to the pseudo-Dirac nature of neutrinos, i.e. v e L and v£ oscillate 
into and v R , respectively, with 45° mixing, depleting the active (to detection) neutrinos 
by one half. 

Since P m depends on the energy E and the SAGE, Homestake and Kamioka experi- 
ments have different energy thresholds, P m can be different for S, H, and /C. Using the 
expression for P m , we have plotted the Am^ - sm 2 (29 efl )/ cos(26> e/ J diagram in Fig. 2 for 
the pseudo-Dirac MSW and the usual MSW effects based on iS exp , H exp and /C exp . We 
have taken (E) = 2.0 MeV, 7.5 MeV and 10 MeV for S, H and /C, respectively. The 
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region which satisfies iS exp , H exp and /C exp within 2a in the standard two generation MSW 
effect (with Eq . fl2j]) ) is shown as the area inside the dotted lines in Fig. 2. The region 
which satisfies iS exp , ^ex P and /C exp within 2a in the pseudo-Dirac case (with Eq.([L9|)) is 
shown as the shaded area inside the solid lines in Fig.2. It is important to emphasize here 
that the two allowed regions in the Am^ - sin 2 (26^)/ cos(20 eM ) plot based on the same 
data, <S exp , H exp and /C exp , do not overlap, even within 2a errors, i.e. the two generation 
pseudo-Dirac neutrinos produce different allowed regions from those based on the usual 
two generation model. As variations of the above pseudo-Dirac neutrino scenario, we 
consider the following two cases: (I) The electron neutrino is a pseudo-Dirac neutrino but 
the muon neutrino is an ordinary Majorana or Dirac neutrino; (II) The electron neutrino 
is a Majorana or Dirac neutrino but the muon neutrino is a pseudo-Dirac neutrino. In 
case (I) there is one MSW resonance region and the neutrino fluxes at the earth depend 
on the corresponding Landau- Zener transition probability P m . By assuming a mixing 
between and z/f with mixing angle 6 e/1 , we have 



vf : i/f 



1 



P m cos 2 6 efM + (1 - P m ) sin 2 9 eix : P m sin 2 ^ + (1 -P m )cos 2 # 



'en 



(23) 



leading to 



s = n 



^-^(l-2P fll )cos(20 e/1 ) 



6 3 



^-^(l-2P m )cos(2^) 



(24) 

When we neglect the energy dependence of the detection rates, there is a region which 
satisfies the three detection rates within 2a uncertainties 



0.40 < 



: (1 -2P m )cos(2fl 



< 0.70. 



(25) 



The allowed region is similar to the case in which both neutrinos are pseudo-Dirac (see 
Eq. (P0"D ) since in both cases there is only one adiabatic MSW transition region. 

In case (II) there are two MSW resonance regions, but in the region P2 the transition 
is extremely non-adiabatic, i.e. P R2 ~ 1. The neutrino fluxes at the earth depend on the 
Landau-Zener transition probability P m . By assuming a mixing between z/f and v e with 
mixing angle 8 efl , we have 



P m cos 2 6 efi + (1 - P R1 ) sin 2 9 etl 



P R1 sin 2 e/U + (1 - P m ) cos 2 9 efl 



(26) 



leading to 



S = H= ~^(l-2P m )cos(2^) 



K 



1 11 
12 + 12 



1 1 



1 - 2P m ) cos(20, 



(27) 



S 



When we neglect the energy dependence of the detection rates, we find a region which 
satisfies the three detection rates within 2a uncertainty as 



0.24 < 



i-^(l-2P fll )cos(20 eM ) 



< 0.35. (28) 



This allowed region is similar to the case of the ordinary MSW result (see Eq . (g2j) ) . Note 
that this case is different from the previous one and the case of both neutrinos being 
pseudo-Dirac because the electron neutrino flux is not depleted in half. 

In summary, the allowed regions in the mass-mixing angle parameter space which 
are obtained from the SAGE, Homestake and Kamioka experiments are very different 
depending on the particle content of the neutrino sector. For example, the region allowed 
when both u e and are pseudo-Dirac neutrinos is very different from the region allowed 
by the analysis based on the usual two generation MSW effect. Consequently, in order 
to pine down the values of Am^ and 6 efl from future solar neutrino experiments, it is 
necessary to have a complete understanding of the neutrino sector, in particular whether 
sterile neutrinos actually exist or not; if they do, what their nature would be and so on. 
Finally, we conclude with short comments on the apparent atmospheric neutrino puzzle 
and the neutrinoless double beta decay. First, the atmospheric neutrino puzzle is that 
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<T>° bs / <|> cal 
<T>obs / ^cal 



0.65 ± 0.08 ± 0.06 Kamioka 

(29) 

0.64 ±0.09 ±0.12 1MB 



where $° bs „ and $f, al „ are the observed and calculated fluxes of atmospheric neutrinos, 
respectively. This puzzle can easily be solved in the scenario in which both v e and are 
pseudo-Dirac, as mentioned in Ref. 0, or in the scenario (II) discussed above as long as the 
constraints Am^ < 10~ 7 eV 2 and 10 _4 eV 2 ^5 Am^ < 10~ 2 eV 2 are met. The upper limits 
are both due to the cosmological argument [|l(J and the lower limit for Am 2 is necessary 
in order to have an oscillation length much shorter than the radius of the earth. In these 
scenarios, is depleted in half simply because of its pseudo-Dirac nature. Furthermore, 
a value Am 2 ~ 10~ 4 eV 2 , which correspond to an oscillation length equal to the earth 
diameter for E ~ 500 MeV, could explain the observed suppression of the flux of low- 
energy muon neutrinos and a value between 0.5 and 1 for the ratio given in Eq.(^). Note 
that, in the case in which both neutrinos are pseudo-Dirac, v e is not depleted because 
the oscillation length is much longer than the radius of the earth due to the cosmological 
limit mentioned above. Secondly, as already discussed in the past |LlJ , neutrinoless double 
beta decay rates are naturally suppressed because they become proportional to m e L 2 , for 
the pseudo-Dirac neutrinos that we have discussed. This implies that non-observation 



9 



of neutrinoless double beta decay cannot automatically lead to the conclusion that the 
electron neutrino is a Dirac particle. 
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Figure Captions 



Fig.l Effective masses squared in matter for the energy eigenstates vi, 1^2, z/3 and 1/4 as 
functions of the matter density p. 

Fig.2 Am^-sin 2 (2# e/i )/cos(2# e/i ) plots. The area inside the dotted line is the allowed 
region in the standard MSW effect, whereas the shaded area is the allowed region 
in the pseudo-Dirac case. The transition at R2 becomes non-adiabatic in the region 
below the solid line in the upper right corner. 
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